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This article investigates the full Boltzmann equation up to second order in the cosmological per- 
turbations. Describing the distribution of polarized radiation by using a tensor valued distribution 
function, the second order Boltzmann equation, including polarization, is derived without relying 
on the Stokes parameters. 
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Introduction 



The properties of the cosmic microwave background 
(CMB) temperature fluctuations depend on both the ini- 
tial conditions set at the end of the primordial inflation- 
ary era, and on their post-inflationary evolution. We 
describe radiation using the kinetic theory, in which its 
properties are encoded in a distribution function. The 
theory of cosmological perturbations around a maximally 
symmetric space-time enables to solve order by order the 
evolution equation for the CMB anisotropics. The lin- 
ear perturbations fail to capture the intrinsic non-linear 
features of General Relativity which can enter both the 
initial conditions and the evolution, and we thus need at 
least the second order if we want to estimate the bispec- 
trum in the CMB and better reveal the physics of the 
primordial universe. Indeed, the bispectrum can only be 
generated either from non-Gaussian initial conditions set 
by inflation or by secondary effects during the subsequent 
evolution, and though standard scalar field slow-roll infla- 
tion predicts negligible amounts of non-Gaussianity 
the latest observations of the CMB 0, 0] indicate that 
there might be a non-vanishing bispectrum. We thus 
need to extend the program followed at first order in per- 
turbations up to second order, that is to build a full set of 
second order gauge invariant perturbation variables and 
derive the perturbed Boltzmann and Einstein equations 
which determine their dynamics. In the fluid limit, the 
gauge issue was studied in Ref. M and gauge invariant 
variables were built in Refs. [H, Tg] up to second order 
in perturbations. This fluid approximation has already 
been used to understand the general form of the bispec- 
trum on small scales generated by evolutionary effects 
in Ref. 0]. As for the kinetic theory, the gauge issue 
was studied in our previous paper ||, and the evolution 
equations through free-streaming and Compton collision 
with free electrons were derived up to second order in 



Refs. 
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Il2j . but it is so far restricted to un- 
polarizcd radiation, which is inconsistent since Compton 
scattering does generate polarization. This letter sum- 
marizes the full second order derivation of the radiation 
transfer, including polarization, detailed in Ref. flij ]. 



Kinetic theory with polarization 

Describing polarized radiation The momentum of pho- 
tons is usually decomposed on an orthonormal basis, that 
is using a tetrad field defined by e Q .ef, = ej^e^g^ = 
r] abl e a .e b = e a ^e h v g^ v = rj ab . We use Greek indices 
(/i, v, p . . . ) for abstract indices and the beginning of the 
Latin alphabet (a, b, c ■ ■ ■ = . . . 3) for tetrad labels. We 
also use i,j,k--- = 1...3 for the spatial type tetrads 
and we reserve the label o for the timelike vector. As 
for the labels associated with the coordinates, we use 
A,B,C,--- = ... 3, and I,J,K,- - = 1 ... 3, with the 
time component labeled by O. 

A momentum can then be written p = p a e a = p°e a + 
p l Gi, and decomposed between the energy p° and the 
spacelike direction unit vector n according to p' 1 = 
P°i e o + n>1 ) with n^eg = 0. This decomposition can 
be used to define the screen projector 



<V/(p) = 9fiv 



e° e° 



(1) 



which projects on the space both orthogonal to e D and 
to the direction n. The radiation is represented by a 
Hermitian tensor valued distribution function also called 
polarization tensor [3j EH satisfying p^F^ v {p a ) = 0. 
The screen-projected distribution function, f^ v {p a ) — 
S?S° F pa (p a ), has four degrees of freedom which can be 
split according to 



(2) 



where e pva = e^e ptllV( j, with e piLVa the space-time fully an- 
tisymmetric tensor, and where we have omitted for sim- 
plicity of notation the dependence in p a . P pVl which en- 
codes the degree of linear polarization, is real, symmetric 
and trace free, as well as orthogonal to e Q and n. It has 
two degrees of freedom, usually described by the Stokes 
functions Q and U. I and V are respectively the inten- 
sity (or distribution function) for both polarizations and 
the circular polarization. The brightness is then defined 
as 



4-7T 



I(p°,n l )( P °ydp° 



(3) 
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and we can define similarly V and V^ v . 
The remaining dependence in n, can be further ex- 
panded in multipoles using projected symmetric trace- 
free (PSTF) tensors, where projected means that they 
are orthogonal to e Q . For instance, using the notation 



and ti 1 = n^e 1 , X can be expanded in 



(4) 



where the Xi t = i t are PSTF, and where convention- 
ally, the lowest multipole, i.e. corresponding to i = 0, 
is noted 2g. A similar expansion can be performed on 
V. As for V a b = T^iiv&a^b ', its non vanishing components 
can be expanded in electric and magnetic type multipoles 
according to 14, 15| 



Pi 
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ijki-2 



TT 



(5) 

where the notation TT denotes the transverse (to n) 
symmetric trace- free part. We have used (..) for the 
symmctrization of indices, and we will also use [..] for 
the antisymmetrization. 

Transformation rules under a change of frame This 
description depends on the tetrad chosen to define the 
polarization tensor, the energy and the momentum direc- 
tion. Under a local boost parameterized by v l = v^e 1 ^, 
the tetrad transforms to e a — e{A 6 a , where the com- 
ponents of the transformation are given by A° = 7 = 
(1 - v.v)-i, A° = -7«j and A*. = 5) + [ 7 2 /(7 + 
Thus, the magnitude and the direction unit vector of the 
photon momentum transform to 



p.e" 



IP 



n.v) 



(6) 



n.v) [ (1+7) 



n.v v 
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This implies that 
c _ c p q & q 

Since the screen-projected polarization tensor transforms 
to 



= n^&LSi . (7) 



Uu(p ,h a )=S ll a S 1 ff a ^p ,n J , 



(8) 



then it can be checked from the decomposition ([2]) and 
the transformation rules that / and V transform 
as scalars, that is 7(p 5 ,n a ) = I{p°,n a ) 1 V{p 5 ,n h ) = 
V(p°,n a ), and that P M „ transforms as f^ y . The differ- 
ential solid angle of the momentum direction transforms 
according to dfl = [7(1 — v.n)] 2 dfi, and thus the trans- 
formation rules of the energy integrated multipoles of the 
brightness can be deduced to be given by 



If 00 
— / dtl [7(1 - v.n)] 2 £ Ij^n^, 

1 J r=o 



(9) 



where A £ = (4tt£!)/[(2£ + 1)!!] and (...) is the notation 
for the symmetric trace-free part. The transformation 
rules of the electric and magnetic multipoles are detailed 

Hi EE I . 
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The Boltzmann equation The evolution of the polariza- 
tion tensor is dictated by the Boltzmann equation 



L[fab] 



5C Q d 
n Oh 



P h V h f, 



cd 



df cd dp h 



dp h ds 



Cab 



(10) 



where L[] is the Liouville operator and C a b the collision 
tensor. It can be shown that the Liouville operator pre- 
serves the decomposition of f^ v in an antisymmetric part 
(V), a trace (/) and a symmetric traceless part (P^v), 
that is 

L[fab] = \L[I]S ab + L[P ab ] + l -L[V]n c e cab . (11) 



Cosmological perturbations 

Perturbation of the metric We assume that, at low- 
est order, the universe is well described by a Friedmann- 
Lcmaitre space-time (FL) with Euclidian spatial sections. 
The most general form of the metric for an almost FL 
universe is 

ds 2 = a(i]) 2 [-(1 + 2$)d?7 2 + 2w/dx J d7 ? + /i/./dx'dar 7 ] , 

where 77 is the conformal time for which the corre- 
sponding index is O, and a(n) is the scale factor. Wc 
define the Hubble parameter by H = a' /a where a 
prime denotes a derivative w.r.t. 77. We work in the 
Poisson gauge though, as already mentioned, a complete 
treatment of perturbations can be performed by defining 
gauge-invariant variables 0, Q, and we perform a 
scalar- vector-tensor decomposition in loj = diB + Bj, 
hu = (1 - 2*)<5/j + 2H U , where 5/ and H u are 
transverse {d I Bi = d 1 Hu = 0), and Hu is traceless 
{H I I = 0). Each of these variables can be split in 
first and second-order parts as W = + ^W^- 2 \ 

and if a quantity has also a background value it is 
noted W. First-order variables are solutions of first- 
order equations, whereas second-order equations will 
involve purely second-order terms, e.g. MA 2 ' and terms 
quadratic in the first-order variables, e.g. [W^ 1 )] 2 for 
which we will omit the order superscript. We neglect 



the first order vector modes (B 



(i) 



and the first order tensor modes {H^j = 0) which arc 
expected to be very small in standard models of inflation. 

Perturbations of tetrads and Ricci rotation coefficients 
The perturbed tetrad can be expressed in function of the 



0) which decay 



background tetrad in the generic form by e a 



(n) 



i?i n)h e 6 , 



e b («) = e a S a ■ The symmetric part R^L is constrained 
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by the normalization conditions of the tetrad and is thus 
expressed in function of the metric perturbations. We are 
free to choose the antisymmetric part iir w , since it corre- 
sponds to the Lorentz transformation freedom (boost and 
rotation) . We require e° ~ dr] [lj| , which is equivalent to 
choose R^g = S 1 -"^ = for any n. We also fix the rota- 



tion by requiring = = 0. This procedure means 
that our tetrad is adapted to observers whose velocity is 
always orthogonal to constant time hypersurfaces. We 
then choose a background tetrad adapted to our coor- 
dinates, i.e. e b c = 5^ /a and e b c = a8 b c , which finally 
enables us to express the perturbations of the tetrad in 
function of the metric perturbations. The perturbations 
of the Ricci rotation coefficients ui a b c = Vbd^eJ^V ^e^ , 
needed to express the covariant derivative in the tetrad 
basis, follows directly from the perturbation of the tetrad, 
and the components used in this paper are = 0, and 



(n) 



,(2) 



ioj 



—d [K ^5 n ., 



(12) 



-<9/$ (2) + 2(2$ - , 
d {I B^-Hf)' + (h$W+*W)5ij 
-(3W$ + 2*')^/./ + 4** / (5 A '(./<5n 



The perturbed Liouville operator 

In practice, we want to express the Boltzmann equa- 
tion in function of ry, since we want to perform an integra- 
tion on coordinates. We should thus multiply Eq. (fTT)]l by 
ds/dn = 1/p . However, there is no point multiplying 
with the full expression of ds/drj since it would then bring 
metric perturbations in the collision tensor. Instead, we 



a/p°, and we 



multiply only by (ds/d-q)^ = (l/p°y 

will use the notation L#[] = L\\ (ds/drj)^ and a sim- 
ilar definition for the collision term. Since we are also 
exclusively interested in the dynamics of the brightness, 
we focus our attention on the energy integrated Liou- 
ville operator defined similarly to Eq. ([3]), and we 
use a similar definition for Cf b . The perturbation of the 
Liouville operator up to second order is of the form 

C*[X,X<U,XW] = C*[X]+£#W[X,XW] (13) 

+1 {r#(2)[^ ; X (V] + c*W[x, x ( V}} , 

where X stands for either X, V or V a b- A similar de- 
composition is performed for the collision term. Since 
£#( 2 )[] contains the terms linear in purely second order 
variables, it has the same functional form as £* < - 1 -'[], and 
consequently we only need the expressions of [] and 
-t report the Boltzmann equation up to second 
order in perturbations. 



We first need the perturbation of the trajectory at first 



order which is [dx 1 /drj) = (p 1 /p°) = rt l ($ + *). 
Then, from the perturbations of the Ricci rotation coef- 
ficients (fT2|) , we deduce the perturbations of the geodesic 
equation which are given by 



dpa 
ds 



(n) 



Xn)rb _ 



bac 



p c p° = . 



(14) 



This is used to obtain the second order evolution equation 
for the energy and the first order evolution equation for 
the direction (the lensing equation) which read 



% r ■ »° 



Hj 2 j ] 



(i) 



h* (2) ' + 2($ - ^)n i 5/$ + 4**' , (15) 
-S ij dj (* + *). (16) 



dn l \ 
drj J 

This is all what is required to obtain from Eq. (fT0|l the 
second order Liouville operator for X = I, V, V a b 



jr*W[X,xW} = 



dX^ 2 
dt] 



nidjX [2) + mX 



(2) 



(17) 



-46 % X 



(2) 



/(2) 



[ dl Bf 



h\ 2 J i 



+8X (1) (n*'dj$ - -8 5%X [($ - ^)n%^ + 2**'] 
a v(l) 

-2S ij j (V + $)-—— -2$£ # « (18) 
on 1 

where it is implied that for X = V a b, according to 
Eqs. (fT0|) and (| 1 1 [) . this expression also needs to be screen 
projected. The notation <5j has also been introduced in 
terms involving X to remind that, due to the symme- 
tries of the background space-time, there is no circular 
or linear polarization at the background level. 



The collision term 

Electrons without bulk velocity We consider the case 
where the free electrons have only a bulk velocity and no 
thermal dispersion in their velocity distribution. Since 
all electrons have the same velocity, we choose to align 
the first vector of the tetrad e Q with this bulk velocity, 
that is to work in the baryons rest frame. In that case, we 
can neglect the terms arising from the recoil of electrons, 
since they will be of order T r /m e ~ 10 -6 around recom- 
bination, where T r is the temperature of radiation, and 
consequently they do not contribute to the bispectrum. 
The Thomson approximation is thus sufficient and , in 
this case, the collision tensor is known and reads 



c*( P a ) 



HO' 

—S.?S l 3 f kl {p' a )-f l3 {p a ) 



la In I 

(19) 
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where r' = an c o~T, with n c the electrons number density, 
and <7t the Thomson cross section. Explicitly, for the 
energy integrated collision term we obtain 



TT 



Si i 

2 3 



1(J 5 



-V(»«) + ivzn' 



(20) 



We clearly see on this expression that the circular 
polarization is not excited and thus remains null if so 
initially. This explains why we did not bother reporting 
the corresponding part of the Liouvillc operator, and 
from now on we neglect it. Now, if the electrons have a 
thermal velocity, which is in fact the case, the previous 
expressions has corrections of order T c /m c ~ T r /m e 
where T e is the temperature of electrons, but as for 
the recoil term, we can discard them as they will not 
contribute to the bispectrum. 

Electrons with bulk velocity If the distribution of elec- 
trons has a bulk velocity relative to the cosmological 
frame, then we just need to transform the result obtained 
in the baryons rest frame to the cosmological frame, using 
the transformation rules ((6]) ((8]) and ((9]) and remembering 
that the energy integration is done according to Eq. ((3]). 
Keeping only the terms which can contribute eventually 
to the second order expansion, we obtain 



-X(n q 



2 S 13 



la* 



-^2 k n k n.v + -l kl n k n l n.v (21) 
^l kl n k n l j (1 - n.v) - ^l k v k 



(6£ u - l H )n k v l + 4Z n.v - Z v.v + TT^v.n)' 



^£ kl n k n\l + 5v.n) - jB kl e m P l v p n m n" 

O 



Vi j {n")(n.v - 1) + ^B/e^ + - 
z 



TT 



+ (6£i k — T ik ) 
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If 



This expression agrees with the results obtained for the 
EE O El (see details in Ref. [T~ 



unpolarized case 



)■ 



The perturbative expansion of the form (|13p can then 
be easily read. The expression of the collision tensor to- 
gether with the Liouville operator obtained in Eqs. p7|) 
and (fT8|) are our key results. It is then possible to shift 
to Fourier space and expand the Liouville operator and 
the collision tensor cither in PSTF multipoles, or in nor- 
mal modes components fl8j | which are better suited for 



numerical integration. Details of these extractions can 
be found in Ref. [l3j]. Hence it offers, in principle, the 
possibility to integrate numerically the Boltzmann hier- 
archy up to second order, though it is expected to be 
very time consuming, and approximate schemes should 
be developed in order to focus on the dominant effects. 
Computing the effects of non-linear evolution has now 
become extremely important since it has been recently 
claimed that a non-zero non-Gaussianity in the CMB 
had been detected Since the precision of the CMB 
measurements is going to increase significantly in the 
forthcoming missions, disentangling the primordial non- 
Gaussianity, which delivers valuable information on the 
primordial universe, from the subsequent non-linear ef- 
fects is a necessary task in order to constrain better from 
inflation our theories of high energy physics. 
Acknowledgements: I thank J-P Uzan for many useful 
discussions on the kinetic theory. 
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